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Abstract. We present an analysis of the optical data 
of the Hubble Deep Field South. We derive F30QW ab , 
F450Wab, F606Wab and F814WAB number counts for 
galaxies in our catalogue (Volonteri et al. 2000): the 
number-counts relation has an increasing slope up to the 
limits of the survey in all four bands. The slope is steeper 
at shortest wavelengths: we estimated 7f300Wab ~ 0.47 ± 
0.05, 7F450W AB ~ 0.35 ± 0.02, 1f606Wab ~ 0.28 ± 0.01 
and ^FiUWAB ~ 0.28 ± 0.01. The color-magnitude rela- 
tions of galaxies shows an initial blueing trend, which gets 
flat in the faintest magnitude bin and the sample con- 
tains a high fraction of galaxies bluer than local sources, 
with about 50% of sources with F814W.4B > 27 having 
(F450W-F606W)ab bluer than a typical local irregular 
galaxy. Morphological number counts are actually domi- 
nated by late type galaxies, while early type galaxies show 
a decreasing slope at faint magnitudes. Combining this 
information with photometric redshifts, we notice that 
galaxies contributing with a steep slope to the number 
counts have 2^1, suggesting a moderate merging. How- 
ever we emphasize that any cut in apparent magnitude at 
optical wavelengths results in samples biased against ellip- 
tical galaxies, affecting as a consequence the redshift dis- 
tributions and the implications on the evolution of galaxies 
along the Hubble sequence. 

Key words: cosmology: observations - galaxies: evolution 
- galaxies: statistics - galaxies: luminosity function, mass 
function 



1. Introduction 

The Hubble Deep Field South (HDF-S) was observed in 
October 1998 by the Hubble Space Telescope (HST). It is 
the southern counterpart of the Hubble Deep Field North 
(HDF-N) and shares its characteristics of depth and spa- 
tial resolution. The HDF-S is a four arcmin 2 survey lo- 



cated at RA 22h 32m 56s, DEC -60° 33' 02", observed 
during 150 orbits with the Wide Field Planetary Cam- 
era 2 (WFPC2). HDF-S images cover a wavelength range 
from the ultraviolet to near-infrared with 4 broad-band fil- 
ters: F300W, F450W, F606W, F814W. Details about ob- 
servations and data reduction may be found in the Hub- 
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ble Deep Field South web page (bttp://www. stsci.edu/ 
ftp /science /hdfsouth/hdfs.html) . 

The main goal of these surveys, HDF-N and HDF-S, 
is to study faint galaxies, in particular to investigate the 
nature of the faint blue galaxies excess and the evolution 
of galaxies (e.g. Ellis 1997). 

Faint galaxy counts have been of paramount impor- 
tance to show that galaxies do evolve with redshift, even 
if the overall scenario and the physical processes that led 
to evolution are still debated. An excess of observed galax- 
ies with respect to a no-evolution model begins to appear 
at B~ 21 - 22 (e.g. Maddox et al. 1990; Jones et al. 1991; 
Metcalfe et al. 1991) and continues to rise up to B> 28 
(Metcalfe et al. 1995) exceeding by a factor 4 to 10, de- 
pending on the choice of q (0-0.5) and on the LF (see 
e.g. Koo & Kron 1992; Ellis 1997). This excess has actu- 
ally undoubtly proved that galaxies evolve. However, the 
analysis of the optical data have not offered a unequivocal 
interpretation: the steepening of the faint end luminos- 
ity function (LF) of galaxies at high redshifts (Lilly et al. 
1995; Ellis et al. 1996) suggests that either the typical lu- 
minosity of galaxies was brighter by ~ 1 magnitude at 
z ~ 0.7 or that there were 2-3 times more galaxies at this 
epoch. Most likely both luminosity and number density 
evolution are acting: luminosity evolution seems to play 
an important role at least on disk galaxies as shown by 
the results obtained by Schade et al. (1995) and Lilly et 
al. (1998) on the CFRS-LDSS sample. On the other hand 
also mergers of galaxies seem to significantly contribute 
to the evolution of both the luminosity function and lu- 
minosity density of the Universe at least down to z ~ 1 
(Le Fevre et al. 2000). Analyzing high resolution ground 
based imaging of blue galaxies at redshift z < 0.7, Colless 
et al. (1994) actually noticed that most of them have a 
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close companion at a projected distance of about 10 kpc. 
This result suggests the presence of merging at z < 0.5, 
which could increase the star formation, explaining the 
blue colors of the sources. 

Disentangling the effects of luminosity and number 
density evolution is fundamental in understanding the 
nature of field galaxy formation and evolution, allow- 
ing a direct comparison with models. For instance, semi- 
analytical models of galaxy formation in hierarchical clus- 
tering scenarios, assuming that galaxies assemble through 
successive merger events from smaller sub-units, make 
concrete predictions about the evolution of merger rate, 
morphological mix and relative redshift distribution (e.g. 
Baugh et al. 1996; Kauffmann 1996; Kauffmann & Char- 
lot 1998). In this scenario the number density of galaxies 
should increase with redshift, their sizes should decrease 
and the fraction of high redshift galaxies (z > 2) in an 
IR-selected sample should be very low. Moreover the mor- 
phological mix should change in favour of late type and 
irregular galaxies which thus would dominate both galaxy 
counts at faint magnitudes and the high redshift tail of 
the redshift distribution, where a lack of elliptical galaxies 
should be observed. Most of these predictions have been 
the object of detailed analysis of data resulting in con- 
flicting results. For instance, some works claim a deficit of 
ellipticals at z > 1 (Zcpf 1997; Franceschini et al. 1998; 
Barger et al. 1999) and a very low fraction (<~ 5%) of 
z > 2 galaxies in a K< 21 selected sample (Fontana et 
al. 1999). On the contrary other authors find a constant 
comoving density of ellipticals down to z ~ 2 (e.g. Bem'tez 
et al. 1999; Broadhurst & Bowens 2000) and a strong clus- 
tering signal of EROs in K-selected samples (Daddi et al. 
2000). HST Medium Deep Survey (MDS), on the other 
hand, indicates that a large fraction of galaxies at moder- 
ate redshift (z ~ 0.5) are actually intrinsically faint and 
exhibit peculiar morphologies suggestive of merging (Grif- 
fiths et al. 1994; Driver et al. 1995; Glazebrook et al. 1995; 
Neuschacfer et al. 1997). The results derived from the mor- 
phological analysis on the HDF-N data (van den Bergh et 
al. 1996; Abraham et al. 1996; Driver et al. 1998) seem to 
confirm the high fraction of apparent irregular/peculiar 
galaxies with respect to the local universe which may be 
the cause of the faint blue galaxies excess (Ellis 1997). 

Williams et al. (1996) stress also that HDF-N galax- 
ies show a strong bluing color trend: while red galaxies 
are compatible with an elliptical galaxy colors, most of 
blue sources are bluer than an irregular galaxy, and these 
colors might be accounted for only considering luminosity 
evolution. 

Thus it seems that the overall evolutionary scenario is 
still missed and the understanding of the contribution of 
the different populations of galaxies to this scenario rather 
confused. 

In this work we present an analysis of number counts, 
colors and morphological distributions of HDF-S galaxies 
in an attempt to reveal and better constraint the evolving 



Table 1. RMS sky values and 5ct magnitude limit for the 
different bands. 
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zero-point 


RMS 








(ADU/pix)xlO" 5 




F300W 


20.77 


1.674 


28.87 


F450W 


21.94 


2.284 


29.71 


F606W 


23.04 


4.126 


30.16 


F814W 


22.09 


2.960 


29.58 



populations of galaxies. The plan of the paper is as follows: 
in Section 2 we summarize the steps leading to the cre- 
ation of the catalogue, the computation of optical galaxy 
counts is explained in Section 3, in Section 4 the technique 
used to determine galaxy colors. In Section 5 we discuss 
our results, both counts and colors, while in Section 6 we 
present and discuss morphological number counts and in 
Section 7, by adding the information given by photomet- 
ric redshifts and number counts models, we disentangle 
the contribution of different galaxy populations and dis- 
cuss the interpretation of redshift distributions. Finally, 
in Section 8 we summarize our results and express our 
conclusions. 

2. Image Analysis 

The samples here used are not extracted from the public 
HDF-S catalogs available on the Web. On the contrary 
they have been independently created beginning from the 
optical HST images in order to have a control of the de- 
tection procedure and of magnitude estimate. A detailed 
description of the procedure adopted and the selection cri- 
teria used to obtain the catalogs in the different bands, as 
well as magnitude estimates, can be found in Volonteri et 
al. (2000) . Here we briefly describe the steps related to the 
derivation of number counts. 

The object detection was performed across the area ob- 
served in optical bands (WFPC2 field) using SExtractor 
(Bcrtin and Arnouts 1996). Images were first smoothed 
with a Gaussian function having FWHM equal to the 
one measured on the images (w 0.2 arcsec). A detection 
threshold of 1 sigma per pixel and a minimum detec- 
tion area equal to the seeing disk (« 0.02 arcsec 2 = 13 
pixel) were adopted to peak up objects. This thresh- 
old corresponds to a minimum signal-to-noise ratio of 
S/Nwf = 1-34 and S/Npc = 0.67 for the faintest sources 
detectable on the WF area and on the PC area respec- 
tively. 

In Table 1 we report for each filter the zero-point (AB 
magnitude, Oke 1974), the sky RMS estimated by SEx- 
tractor and the corresponding 5(7 limiting magnitude for 
a point source. 

"Pseudo-total" magnitudes were estimated using the 
method of Djorgovski et al. (1995) and Smail et al. (1995) 
on the basis of the following steps: 
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— the SExtractor isophotal corrected magnitude has been 
assigned to large sources, i.e. those sources having an 
isophotal diameter D iso > 8\ and to those sources 
nagged by SExtractor as "blended" ; 

— the aperture corrected magnitude (estimated within 
9 1 and then corrected to 82, being 8\ < 8 2 ) has been 
assigned to small sources with Di SO < 8\\ 

8\ has been defined as the minimum apparent diameter of 
a galaxy having an effective diameter r e = 10 kpc. Here- 
after we use a A = cosmology, with go = 0.5 and Hq = 50 
km s _1 Mpc _1 unless differently specified. With this choice 
8\ = 1.2 arcsec. 82 is the diameter corresponding to the 
area for which 90% of the bright sources have a smaller 
isophotal area (see Volonteri et al. 2000 for a complete de- 
scription of procedure, tests and results). This procedure 
was applied independently for each band. 

A measure of the detection reliability is necessary in 
order to evaluate the number of spurious sources included 
in the sample. We treated this problem statistically, in the 
hypothesis that noise is symmetrical with respect to the 
mean sky value. Operationally we first created for each fil- 
ter a noise frame by reversing the original images to reveal 
the negative fluctuations and to make negative (i.e. unde- 
tectable) real sources (Saracco et al. 1999). Then we run 
SExtractor with the same detection parameter set used 
to search for sources in the original images detecting, by 
definition, only spurious sources. In Figure |l| the magni- 
tude distribution of spurious sources obtained on the WF 
area and the PC area in the F606W band are shown re- 
spectively. Applying a S/N=5 cut off is sufficient to re- 
duce the spurious contamination to a negligible fraction 
(4%) on the WF area, after removing the edges of the 
images with lower sensitivity. On the contrary such a cut 
off is not able to reduce spurious detections to a reason- 
able level on the PC area, contamination being more than 
35%. Thus due to such a large number of spurious sources 
which would be introduced by the PC data, we restricted 
the selection of sources to the central WF area only cor- 
responding to 4.38 arcmin 2 . On this area 450, 1153, 1694 
and 1416 sources were selected accordingly to the above 
criteria in the F300W, F450W, F606W and F814W band 
respectively. 

3. Optical Galaxy Counts 

To derive galaxy counts, we first removed stars from the 
sample by using the SExtractor star/galaxy classifier. We 
defined as stars those sources brighter than F814W.4B = 
22 and having a value of the SExtractor "stellarity" in- 
dex larger than 0.9. This choice will tend to underesti- 
mate stars both at faint magnitudes where no classifica- 
tion is considered, and at bright magnitudes where some 
fuzzy stars could be misclassified as galaxies. On the other 
hand this will ensure that our galaxy sample is not biased 
against compact galaxies. The star "cleaning" procedure 
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Fig. 1. Spurious sources, defined as detections on the re- 
versed (i.e. multiplied by -1) F606W image. The upper 
panel refers to the WF field, the lower panel to the PC 
field. The solid line represents the initial detections, the 
dashed line represents spurious sources left after applying 
our selection criteria (i.e. a S/N=5 cut off, after removing 
the edges of the images) . 



has classified and removed 14 stars with F814Wab < 22 
in agreement with the number of stars found in the HDF- 
N by Mendez et al. (1998) to this depth and in excess by a 
factor of two with respect to the prediction of the galaxy 
model of Bahcall & Soneira (1981). 

Completeness correction for faint undetected sources 
strongly depends on the source apparent spatial structure 
and on their magnitude. The high resolution of HDF-S im- 
ages allows the detection of sub-galactic structures, such 
as HII regions. Moreover the light observed from a high 
fraction of the galaxies in the HDF-S is emitted in the 
UV and F450W pass-bands, so that the observed mor- 
phology is strongly affected by star formation episodes. 
These features imply that "typical" profiles of galaxies 
are not able to well describe the shapes of a lot of galaxies 
in the HDF-S. Thus, in order to reproduce the manifold of 
shapes which characterizes sources in the HDF-S, follow- 
ing Saracco et al. (1999), we generated a set of simulated 
frames by directly dimming the original frames by vari- 
ous factors while keeping constant the RMS. This proce- 
dure has allowed us to avoid any assumption on the source 
profile providing an artificial fair dimmed sample in a real 
background noise even if it cannot take into account differ- 
ences due to evolution. We thus define the correction fac- 
tor c as the mean number of dimmed galaxies which should 
enter the fainter magnitude bin over the mean number of 
detected ones. In Table. 2 we report the raw counts n r , 
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Table 2. Differential number counts in F300Wab, F450Wab, F606W A b and F814W ab bands over the inner 4.38 
arcmin 2 of the WF area. 
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the completeness correction factor c, the counts per square 
degree corrected for incompleteness N and their errors un ■ 
Colley et al. (1996) suggested that galaxies in the 
HDF-N may suffer from a wrong selection. High redshift 
galaxies on optical images have a lumpy appearance: first 
the redshift moves the ultraviolet rest-frame light into the 
optical, so galaxies are observed in UV rest-frame, 
where star- forming regions are more prominent; second 
the fraction of irregular galaxies is higher than locally 
(van den Bergh et al. 1996, Abraham et al. 1996) and 
a large number of galaxies may display asymmetry and 
multiple structure. We treated this feature analyzing our 
sample in F814W and F450W. The F814W-band cata- 
logue should suffer less from the effects described above, 
being selected in the reddest filter, the vice versa is true for 
the F450W-band catalogue (but see Volonteri et al. 2000 
for a more detailed analysis). About 20-30% of sources 
in F450W-band catalogue have separation < 1 arcsec. 
We therefore analyzed these sources, by cross-correlating 
F814W and F450W catalogues. In the F450W-band cata- 
logue we selected pairs with a separation < 1 arcsec which 
were not included in the F814W-band catalogue. These 
objects were single sources split in the F450W-band (with 
F450W J 4s ~ 27-29), corresponding to a single detection in 



the F814W-band. We then used SExtractor on the F450W 
frame, after choosing a higher deblend_mincont= 0.1. 

87 sources, with 21 <F450W J 4s < 26, corresponding 
to about 7% of the whole sample, were then considered 
as single galaxies. In Figure || counts obtained with the 
uncorrected catalogue are shown as stars, while corrected 
counts are shown as empty circles. In the faintest bins the 
correction is within the error, accounting for less than 10% 
of the counts. 

Errors were obtained by quadratically summing the 
Poissonian contribution a nr = \prC/ of raw counts, the 
contribution due to clustering fluctuation 

a^co^N (1) 

being lo{9) = A^O^^^ 1 ^ the angular correlation function, 
and the uncertainty on the correction factor 

a- c = ^l/fc^-c) 2 (2) 

where k is the number of frames dimmed for each dim- 
ming factor. Assuming that the amplitude evolves 
with magnitude on the basis of the relation logA w = 
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Fig. 2. Comparison of F450W-band number counts un- 
corrected (stars) or corrected (empty circles) for "over- 
splitting" of sources. The faintest bin are unaffected by 
the correction, since it is within the error. 



-0.3 x mag + const (Brainerd et al. 1994; Roche et al. 



1996) we derived 






log 4^300 Was) 


= 4.438 - 


0.3F300W AB 


log A^FAWWab) 


= 4.314- 


0.3F450W A b 


log AJF606W AB ) 


= 4.438 - 


0.3FQ06W A b 


log A^FSUWab) 


= 4.161- 


0.3F8UWab 



The number counts here derived in the F300W and 
F450W bands and in the F606W and F814W bands are 
shown in Figures [| and |] together with those from the lit- 
erature. The relation between counts and magnitude may 
be written as 



d log N 
dm 



Jim + c\ to < ml 
J2 m + C2 to > m! 



with 7i > 72. The knee ml is at B«25 (Lilly et al. 1991, 
Metcalfe et al. 1995), and the value of 71 varies between 
0.4 and 0.6 according to the band and analogously for 72 
it varies between 0.2 and 0.5. We estimated both 71 and 
72 in F450Wab, F606Wab and F814Wab, obtaining for 
7i : 1i,F450W A b ~ °.4 ± -!> lt,Feoew AB ~ 0.34 ± 0.1 and 
1i,F8UW A b ~ 0.62 ±0.1, and for 72: J2,F450W AB ~ 0.19 ± 
0.01, j2,F6oew A B ~ 0.19±0.1 and 72,_F8i4M/ AB ~ 0.19±0.1. 

The most notable feature appears in the F300W-band 
counts: these are described by a slope Jf30ow A b = 0.47 ± 
0.05 which is much steeper than the value "/F30ow AB ~ 
0.15 derived by Williams et al. (1996) and Pozzetti et 
al. (1998) on the HDF-N data in the same magnitude 



Fig. 3. A comparison of number counts computed in this 
work with counts obtained in literature in Vab and Iab- 
In particular Williams et al. (1996, will96) refer to the 
HDF-N, not corrected for incompleteness. Other symbols 
correspond to Arnouts et al. (1999, arn99), Gardner et 
al. (1996, gard96), Smail et al. (1995, smail95), Lilly et 
al. (1991, lilly91). Our best fit slopes are 7_f606Wab = 
0.28 and Jfsi4W ab = 0.28, in agreement with previous 
works, d\ogN/dm w 0.27- 0.32 (e.g. Smail et al.(1995), 
Arnouts et al. (1999)), while Williams et al. (1996) find 
dlogN/dm = 0.17 for F606Wab=26-29, dlogN/dm = 
0.35 for F606Wab=23-26 and d\ogN/dm = 0.18 for 
F814W A i3=26-29, dlogN/dm = 0.31 for F814W A s=23- 
26. 



range. Such a steep slope, which agrees with the find- 
ings of Hogg et al. (ju ~ 0.5; 1997) and Fontana et al. 
(ju ~ 0.49; 1999), does not depend on a possible over- 
estimate of the incompleteness in the faintest magnitude 
bins: the same slope is actually described by counts at 
magnitudes U300 < 26 where the sample is 100% complete. 
Moreover the F300W-band counts do not show evidence 
of any turnover or flattening down to F300Wyi_B=27, con- 
trary to what is claimed by Pozzetti et al. (1998). We will 
discuss this issue in Section 5.1. The amplitude and slope 
of the number counts in the whole range of magnitudes 
in the other optical bands are in a good agreement with 
those previously derived by other authors. We estimated 
r 7F450W AB ~ 0.35 ± 0.02, 7f606Wab ~ 0.28 ± 0.01 and 
1F8UW AB ^ 0.28 ±0.01 for the F450Wab, F606Wab and 
F814Wa_b counts respectively. 

4. Optical Galaxy Colors 

In order to measure colors unbiased with respect to the se- 
lection band, we created a metaimage by summing all four 
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Fig. 4. As in Figure ||, in Bab and Uab bands. The com- 
parison with previous works includes Williams et al. on the 
HDF-N (1996, will96), Arnouts et al. (1999, arn99), Gard- 
ner et al. (1996, gard96), Metcalfe et al. (1991, met91), 
Tyson et al. (1988, ty88), Jones et al. (1991, jon91). 
Slopes are jp 45qwab =0-35 e jjj = 0.47. Tyson et al. 
(1988) find dlogN/dm = 0.45, Arnouts et al. (1999) 
dlogN/dm = 0.31 in the F450W A b- In F300Wab our 
result is consistent with dlogN/dm w 0.4 — 0.6 (Hogg 
et al. (1997), Koo (1986)). Williams et al. (1996), find 
dlogN/dm = 0.16 for F450WAB =26-29, dlogN/dm = 
0.39 for F450Was=23-26 e dlogN/dm = 0.05 for 
F300Was=26-28, dlogN/dm = 0.40 for F300Was=23- 
26, Pozzetti et al. (1998) find dlogN/dm = 0.135. 

frames, after normalizing them to the same rms sky noise 
(see also Volonteri et al. 2000). We run SExtractor in the 
so-called double image mode: detection and isophote 
boundaries were measured on the combined image, while 
isophotal magnitudes were measured on F300W, F450W, 
F606W, F814W images individually. Using this pro- 
cedure (Moustakas et al. 1997) both object de- 
tection and isophote determination are based on the 
summed image, and isophotes are not biased towards any 
of the bands. 

Then we cross-correlated the catalogue obtained from 
the combined image with F814W and F300W samples se- 
lected according to our criteria (see Section 2): we ended 
up by having two sample selections. Starting from each 
sample, we assigned a lower limit in magnitude to sources 
undetected in any of other bands. The limiting magnitude 
is the 5cr isophotal magnitude within the isophote mea- 
sured in the combined image. In Figures [|, ^, fj| ^ the 
color-magnitude diagrams of the two samples are shown 
together with the median locus. The error bars are the 
standard deviation from the mean of the values in each 



Fig. 9. Fraction of sources in the F814W selected sam- 
ple with (F450W— F606W),4s bluer than a typical irreg- 
ular galaxy (Coleman, Wu, Weedman, 1980). The solid 
line refers to sources detected both in F814W and F450W 
bands, the dashed line histogram includes also lower lim- 
its to magnitudes for sources undetected in the F450W 
band. The dot-dashed histogram includes also galaxies un- 
detected both in F606W and F450W bands, with lower 
limits in both magnitudes. 

bin. The most notable feature is the initial trend towards 
blue colors followed by a flattening in the last bins. 

We then compared our colors with colors measured 
from spectra by Coleman, Wu & Weedman (1980, CWW 
hereafter) after convolution with throughput curves for 
the WFPC2 filters. In each magnitude bin, we estimated 
the fraction of sources with (F450W— F606W)ab bluer 
than an irregular galaxy at z = (i.e. F450Wab— 
F606Wab <0.18). As shown in Figure |, at F8UW A b > 
27 about 50% of sources have (F450W-F606W)ab bluer 
than a typical local irregular. For the F300W-selected cat- 
alogue we also estimated the percentage of galaxies with 
(F300W— F450W)^b bluer than a typical irregular: at 
F300W^b=26 this percentage is almost 80%. 

We will discuss these features in Section 5.2. 

5. Discussion 

5.1. Number Counts 

Our results show a decreasing slope at redder wavelengths 
for the whole sample, in good agreement with other works. 
The very faint end is flatter than 0.2 in F450Wab, 
F606Wab, F814Wab, while in F300WAB the slope is 
steeper. The difference with Pozzetti et al. F300W.4B 
counts is likely due to a different selection of the cat- 
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Fig. 5. Left panel: (F450W— F814W)ab vs F814W A b for the F814W-selected sample: in every magnitude bin the 
median colour is shown, error bars are the standard deviation of the colour distribution. Asterisks represent sources 
detected both in F450W and F814W bands, while circles mark sources undetected in the F450W band: their colour 
is a lower limit, defining the F450W A s magnitude as the 5<r flux within the isophote on the combined UBVI image. 
Right panel: (F450W— F606W) J 4s vs F814WAB for the F814W-selected sample: asterisks represent sources detected 
in all F450W, F606W and F814W bands, circles mark sources undetected in the F450W band and triangles sources 
undetected both in F450W and F606W bands (we assigned to those sources an arbitrary value (F450W— F606W) J 4s = 
-1.25). 



alogue: they selected the sample in a red band 
(on a summed F606W+F814W image) to derive F300W- 
band counts and used the isophotal magnitude instead of 
a pseudo-total magnitude. They assume that galax- 
ies fainter than F606W.4B < 26 have approximately 
a (F300W-F606W)ab color not bluer than the brighter 
galaxies, not biasing the red-selected sample against very 
blue objects. 

In Figure we show the (F300W-F606W) AB color as 
a function of F300W ab isophotal magnitude (as Pozzetti 
et al.) for our F606W-selected sample, within the limit- 
ing F606Wab (i.e. « 30.2) and for the F300W-selected 
one (Figure |ll|) . The colors for the F606W-selected sam- 
ple are obtained as described in Section 4. In order to 
obtain a sample with characteristics similar to Pozzetti 
et al. sample, we considered for F300Wab magnitude the 
flux within the isophote in the metaimage. 

The solid oblique line represents the (F300W 
-F606W)ab color of a F606W A s = 29 galaxy (the 
F606W-band limiting magnitude of Pozzetti et al.): all 
galaxies in the hatched area would be not counted since 
they are fainter than F606W A b > 29, the vertical line is 
the F300W j4 b=29.5 limit. The F606W A s limiting magni- 
tude biases the sample against blue objects, that is the 
last bins are affected by censored data (see next Sec- 



tion), thus making the (F300W— F606W)ab color ap- 
pear redder. Considering the bluing trend observed un- 
til F300W AS =29, and the fact that the F606W AB lim- 
iting magnitude biases the sample against blue objects, 
this implies that F300W-band counts recovered from a 
red-selected sample are affected by this color bias and con- 
sequently by incompleteness in the last bins. The flatten- 
ing found by other Authors may be mainly caused by the 
influence of such incompleteness (noise being equally im- 
portant for faint galaxies, be them at low or high rcdshift) 
and not to the "crossing" of the Lyman break. 

5.2. Colors 

As explained in Section 4 we selected two samples: the 
F814W-selected one should trace the global features of the 
HDF-S sample and the F300W-selected sample evidencing 
star-forming galaxies at moderate z. 

The F814W-selected catalogue has by far the biggest 
number of sources, so that the number of galaxies with 
lower limits in the other bands is rather high. We applied 
a linear fit to the color-magnitude relation by making use 
of "survival analysis" (Avni et al. 1980, Isobe et al. 1986), 
in order to avoid uncertainties introduced by the heavy 
influence of censored data. We applied a linear regres- 
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Fig. 6. Left panel: (F606W— F814W)ab vs F814W ab for the F814W-selected sample: asterisks represent sources 
detected in F606W and F814W bands, circles mark sources undetected in the F606W band. Right panel: 
(F300W-F450W)ab vs F8UW A b for the F814W-selected sample: asterisks represent sources detected in all F300W, 
F450W and F814W bands, circles mark sources undetected in the F300W band or F450W and triangles sources 
undetected both in F450W and F300W bands (we assigned to those sources an arbitrary (F300W— F450W)ab = 2.5). 
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Fig. 7. Left panel: (F450W-F814W)ab vs F300Wab for the F300W-selected sample: in every magnitude bin the 
median colour is shown, error bars are the standard deviation of the colour distribution. Asterisks represent sources 
detected in F450W and F814W and F300W bands, while circles mark sources undetected in the F450W or F814W 
band. Right panel: (F450W— F606W)ab vs F300Wab for the F300W-selected sample. Asterisks represent sources 
detected in F450W and F606W and F300W bands, while circles mark sources undetected in the F450W or F606W 
band. 



M. Volonteri: Interpreting the optical data of the Hubble Deep Field South 



9 




21 22 23 24 25 26 27 28 29 
F300W >n 



X3 

o 

o 

CJ 



& 

o 
o 

CO 
[x, 



o 

LO 
"1- 

* 

o 
o 

CO 

El. 



21 



$Wi I. 




n ^OD O oO °n 

* o o 



22 



23 



24 25 26 
F300W >o 



27 



28 



29 



Fig. 8. Left panel: (F606W-F814W)ab vs F300Wab for the F300W-selected sample. Asterisks represent sources 
detected in F606W and F814W and F300W bands, while circles mark sources undetected in the F606W or F814W 
band. Right panel: (F300W— F450W)ab vs F300Wab for the F300W-selected sample. Asterisks represent sources 
detected in F450W and F300W bands, while circles mark sources undetected in the F450W band. 



sion based on Kaplan-Meier residuals. This analysis has 
suggested an almost flat relation for (F450W— F814W)ab 
vs F814Wab, (F450W-F606W) j4 b vs F814Wab and 
(F606W-F814W)ab vs F814Wab for F814WAB > 24, 
while (F300W-F450W) j4 b vs F814W A b did not reach 
convergence. Smail et al. (1995) noticed a similar tendency 
on their sample limited at R=27. After an initial bluing 
the median V— R color gets redder, V— I gets flat, while 
R— I keeps on following a bluing trend. 

This trend may suggest the presence of a flat spec- 
trum population, whose colors seem to saturate but whose 
contribution is more and more important at faint mag- 
nitudes, as confirmed by the rising fraction of galaxies 
with (F450W— F606W)ab bluer than a typical irregular 
galaxy. Broadly speaking galaxy colors are dominated by 
blue sources, as noted by Williams et al. (1996)in the 
HDF-N, that is the median color is always bluer than 
typical local samples. The F300W-selected sample shows 
a very blue (F300W-F450W) as color, though the rela- 
tion (F300W-F450W)ab vs F300Wab seems almost flat. 
These two features suggest a considerable contribution 
by flat spectrum sources. Comparing the median (F300W 
-F450W)ab color with Bruzual k Chariot (1993) and 
CWW predicted colors, the main contribution may be 
due to sources with z > 0.5, while (F450W-F814W)ab, 
(F606W-F814W)ab and (F450W-F606W)ab are com- 
patible with those of local irregular galaxies. 



6. Morphological Number Counts 

We performed a morphological analysis of the F814W- 
selected sample, limited to F814Wab=25 (Volonteri et al. 
2000). We chose a quantitative classification of galaxies 
morphology, following Abraham et al. (1994, 1996). For 
each galaxy we measured an asymmetry index (A) and a 
concentration index (C a br)- The former is determined by 
rotating the galaxy by 180° and subtracting the resulting 
image from the original one. The asymmetry index is given 
by the sum of absolute values of the pixels in the residual 
image, normalized by the sum of the absolute value of the 
pixels in the original image and corrected for the intrinsic 
asymmetry of the background. The concentration index is 
given by the ratio of fluxes in two isophotes, based on the 
analysis of light profiles. 

We determined number counts by splitting the sam- 
ple in three bins: E/S0, spirals and irregular/peculiar/ 
interacting according to our quantitative classification. In 
every bin the error is estimated as explained in Section 
3, considering only contributions due to Poissonian and 
clustering fluctuation, since according to our simulations 
the sample at F814Wab=25 is complete. 

In Figure |l2| the decrease in differential counts for 
early-type galaxies at F814Wab >22.5, and the steeper 
slope for spiral galaxies and irregular/peculiar/interacting 
ones is clear. They are described by a slope ^i rr = 0.43 ± 
0.05, ^spiral = 0.37 ± 0.05, we discuss this feature further 
in Section 7. 
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Fig. 10. (F300W-F606W)ab color as a function of 
F300Wab isophotal magnitude for the F606W-selected 
sample. Circles represent galaxies with a lower limit 
in the F300W^b magnitude and consequently in the 
(F300W-F606W)ab color. The solid lines represent the 
(F300W-F606W)ab color of a F606W AB =29 galaxy 
(F606W-band limiting magnitude considered by Pozzetti 
et al. 1998) and the F300Wab=29.5 limit. The apparent 
trend towards redder colors for F300Wab >28 is due to 
the F606W-band limiting magnitude. 
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Fig. 11. (F300W-F606W)ab color as a function of 
F300WAS isophotal magnitude for the F300W-selected 
sample. The solid line represent the (F300W— F606W),4s 
color of a F606W j 4b=29 galaxy (F606W-band limiting 
magnitude considered by Pozzetti et al. 1998). 




F814W AB 

Fig. 12. Morphological number counts: early type galax- 
ies are shown as circles, spiral galaxies are shown as pen- 
tagons and peculiar galaxies as squares, the total counts 
are shown as filled circles. Late type galaxies have steep 
slopes, jirr = 0.43 ± 0.05, ^spiral = 0.37 ± 0.05. Lines 
show number counts model for a flat cosmological model 
(q = 0.5, H = 50 km s^ 1 Mpc~\ A = 0) and the LF of 
Marzke et al. (1998): the models by types are obtained by 
keeping the parameters which best fit the total counts. 

Abraham et al. (1996) and Driver et al. (1998) stud- 
ied the morphological number counts for all galaxies with 
F814Wab <26.0 detected in the HDF-N. Our results are 
in agreement with theirs: number counts are dominated 
by late type galaxies and early type sources show a flat 
curve. Results for late type galaxies are compatible with a 
strong evolution in number or with a non negligible frac- 
tion of sources being at moderate redshift, as suggested 
also by the analysis of colors. This result fits in Marzke 
et al. (1997) findings about the LF of galaxies as a func- 
tion of their color: the slope of the LF faint end is steeper 
at low luminosities for the bluest galaxies (Mg w —19.4, 
a < —1.7 for galaxies bluer than B— V=0.4). Actually by 
integrating this LF, within z = 0.5, the number density 
of galaxies bluer than an Irregular galaxy is compatible 
with the number found in our sample. The flattening of 
the curve for early type sources may be due to a decrease 
in density at higher redshift of elliptical galaxies, but it is 
also expected by the a > — 1 slope of the LF of early type 
galaxies (e.g. Marzke et al. 1998). Number counts models 
nevertheless show that an a > —1 LF does not account 
for such a steep decrease (Figure |l^) . 

Brinchmann et al. (1998), Im et al. (1999) and Driver 
et al. (1998) studied the redshift distribution for different 
morphological types in deep surveys. Brinchmann et al. 
(1998) analyzed 341 galaxies at z < 1 from the Canada- 
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France Redshift Survey (CFRS, Lilly et al. 1995, limited 
at F814Was=22.5), and Autofib/Low-Dispersion Survey 
Spectrograph (LDSS, Ellis et al. 1996, limited at b 3 = 24). 
They find the distribution of elliptical galaxies peaked at 
z w 0.5, while spiral galaxies have a shallow distribution 
and the excess of peculiar sources begins at z « 0.4. Ac- 
cording to their conclusions, the number of regular galax- 
ies is compatible with a passive evolution model, while 
the excess of peculiar ones suggests an active luminosity 
evolution and-or number evolution. Im et al. (1999) for a 
sample of 464 galaxies limited at 1= 21.5 do agree with 
Brinchmann et al. (1998). Driver et al. (1998) analyzed 
HDF-N sources, with photometric redshifts, and underline 
the excess both in spiral galaxies at z = 1.5 and irregu- 
lars at z > 1, suggesting number evolution and passive 
evolution for dwarf LSBG at z < 1 giving the density of 
irregular galaxies at low redshift. 



2.5 




7. Constraining the Redshift of Galaxies from 
Photometric Data 

In a deep survey, such as the HDF-S, the contribution 
of low redshift foreground sources should be disentangled 
from the contribution of high redshift ones. We first tried 
to split the two populations with very simple considera- 
tions, regarding galaxy colors and counts. 

We compared colors of HDF-S galaxies with those pre- 
dicted by CWW spectra and by synthetic spectra of 
Bruzual & Chariot (1993, BC hereafter). While BC ac- 
count for evolution, the CWW spectra take into ac- 
count K — corrections. In Figures [T^-^4| apparent colors 
from CWW and BC are shown. We limited our analysis 
at z < 2 to avoid relying too much on the UV part of 
spectral energy distributions which has a strong influence 
at higher redshift. 

By comparison with colors predicted by CWW or BC 
spectra, the observed very blue galaxies (med(F450W J 4s — 
F606WAS) = 0.06, med(F606WAB- F814Wab) = 0.24) 
are compatible with two very different ranges in redshift: 
z < 0.2 or z > 1.5 (and perhaps beyond our analy- 
sis limit). The slope of morphological number counts is 
in good agreement with this hypothesis, assuming that 
most of the faint late type galaxies are at moderate red- 
shift. In the hypothesis of the moderate redshift popula- 
tion (z < 0.2) we computed number counts as a func- 
tion of optical colors, which should reflect rest-frame col- 
ors for this range in z. We divided the F814W-selected 
sample in a "red" subsample, composed by sources red- 
der than the median colour of the sample (F450Wab— 
F606WyiB=0.45) and a "blue" subsample of the remaining 
sources. We limited our analysis at F814W J 4B=26, where 
the catalogue is complete, according to our simulations. 

The curve for blue sources is much steeper than the 
red one, as shown in Figure [l5|. We estimated ^uue ~ 
0.49 ± 0.01. This value is therefore compatible with the 
previous hypothesis, that a non negligible fraction of the 



Fig. 13. (F450W-F606W)ab color from Coleman, Wu 
e Weedman (1980) spectra, convolved with HST filter re- 
sponses. K — corrections have been applied, but no evolu- 
tion has been considered. The solid line refers to irregular 
galaxies, the dot-dashed line to disk galaxies, the dashed 
line to ellipticals. 



blue galaxies may have a moderate redshift, contributing 
with an almost Euclidean slope to the counts or with num- 
ber evolution, with very faint galaxies merging at high red- 
shift. However we cannot rule out the hypothesis that such 
a steep slope can be due to earlier-type evolving galaxies 
which move into the blue sample at high redshift. In such 
a case the steep counts described by the faint blue sam- 
ple would not require a significant rate of merging to be 
justified. 

In order to check the low-redshift solution we com- 
pared the characteristics of those galaxies to the predic- 
tions of a local luminosity function. The brightest of the 
very blue sources would have Mb > —18 if placed at 
z = 0.2 (assuming _ffo=50 km s -1 Mpc -1 and a 
(F450W-F814W)ab=1). 

Considering the local luminosity function (LF) of blue 
galaxies as estimated by Marzke et al. (1997), described 
by a = —1.8 and M B « —19.5 and integrating within 
z = 0.2 and Mb = —18, the number density found 
is compatible with all very blue galaxies (i.e. galax- 
ies with (F450W-F606W)as and (F606W-F814W) j4 s 
bluer than a local Irregular galaxy) being at such a low 
redshift. The first raw technique used to constrain the red- 
shift of sources suggests that a high fraction of faint galaxy 
is composed by low redshift (z < 0.2) sources, nevertheless 
the "high redshift solution" is not ruled out, since steep 
counts may be due also to merging and our color analysis 
is very naive. 
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Fig. 14. As in Figure [T^ but considering Bruzual & Char- 
lot (1993) of synthetic stellar evolution spectra, with a 
Miller- Scalo IMF. The star formation rate is exponentially 
decaying with timescales 1 Gyrs and 5 Gyrs for E and Sb 
respectively, and is constant for Im. 
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Fig. 15. Number counts for blue and red galaxies (i.e. 
(F450W-F606W)ab < 0.45 or vice versa). The blue 
sources have a very steep slope (jbiue ~ 0.49 ±0.01), while 
red sources have a shallower end. 



We then made our analysis better by means of the 
photometric redshift technique. We computed photomet- 
ric redshifts for the galaxies in the F814W catalogue by 
means of the public code hyperz (Bolzonella et al. 2000). 
The method basically consists in a comparison between 



the observed magnitudes and the photometry expected 
from a set of Spectral Energy Distributions (SEDs). The 
efficiency of the method relies on the presence of strong 
spectral features, in particular the 4000 A break and the 
Lyman break for galaxies at high redshift. To reduce the 
uncertainties estimating z p hot, the set of filters must be 
able to identify these characteristics, spanning a wide 
range of wavelength. The photometric redshift z p hot is 
computed by a \ 2 minimization, considering the whole set 
of possibilities with different combinations of the involved 
parameters. 

In this procedure we decided to use the SEDs built 
from the Bruzual & Chariot's synthetic library (GIS- 
SEL98, Bruzual & Chariot 1993), rather than the ob- 
served mean local spectra by CWW, because tests on 
HDF-N indicate a slightly lower dispersion at z < 2. In 
any case, CWW spectra lead to not very different conclu- 
sions. Hence, we selected 5 spectral types, characterized 
by different star formation rates, matching the observed 
colors of the morphological galaxy sequence; a subsample 
of ages allowed by the GISSEL library is also considered. 
Moreover, we took into account the possible presence of 
dust applying the reddening law by Calzetti et al. (2000) 
for different values of Eb-v- The flux decrement produced 
by intervening neutral hydrogen is computed following the 
recipe by Madau (1995). Only spectra with solar metal- 
licity are considered here, because the metallicity can be 
regarded as a secondary parameter. 

The accuracy of photometric redshift estimate ob- 
tained from the considered set of filters can be studied 
by means of simulations on synthetic catalogue. 

Some degeneracy can be found (Bolzonella et al., 
2000), with a non negligible part of galaxies lying at z ~ 1 
- 2.4 incorrectly located at low redshift. These galaxies 
show frequently a probability function with several com- 
parable peaks, at low and high z, due to the degeneracy 
in the parameter space. Near infrared data can in princi- 
ple avoid these uncertainties. Nonetheless, also with the 
available set of filters, we can guess that the objects with 
Zphot ^ 1 really belong to this redshift, because high-z 
objects are rarely misidentified as low- 2 galaxies. 

Photometric redshifts for our catalogue indicate that 
all very blue galaxies are at z p hot > 1 (nied(z p ^ ot ) = 1.6). 
This feature, along with their steep slope (see Section 5), 
suggests a high merging rate for galaxies at z p hot > 1- Vice 
versa very red sources have a shallower redshift distribu- 
tion (med(z p hot)=0-6), though at F814Wab > 25 most 
sources have a high redshift. 

Similarly, at F814W.4B <25, late type galaxies have 
a higher median redshift (med(z p / lo t)=l.l) than early 
type ones (mcd(z p ^ ot )=0.6) and the redshift distribu- 
tion (Figure |l6|) is similar to Driver et al. (1998) one, 
though the number of sources at z p hot > 1.5 is lower. 
There is a point to be emphasized about redshift distri- 
butions and morphological types: the detection of galax- 
ies at high redshift is seriously biased against elliptical 
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Fig. 16. Morphological redshift distributions for 
F814Was <25. 

galaxies because of K — corrections in the F814W-band. 
At z 1, early type galaxies have K — corrections (» 1 
mag) greater than late type K~ corrections , making them 
fainter. It means that a cut in apparent magnitude, such 
as F814Wab <25, biases the sample toward galaxies with 
lower K — corrections. To skip this problem, we tried to 
select a volume-limited sample, based on our photomet- 
ric redshifts. When selecting all galaxies with z p j lo t < 1, 
we found only 15% to be with F814Was <25, that is 
with a reliable morphological classification according to 
our method. Vice versa it means that a cut off in appar- 
ent magnitude does not offer a fair sample for a redshift 
distribution of the different morphological types. 

Probably also the apparent decrease in early-type 
galaxies number counts, as shown in Figure [l^, is due to 
this bias. 

We compared our number counts with different mod- 
els, by using galaxy counts models by Gardner (1998). We 
adopted a flat cosmological model (go = 0.5, H = 50 
km s _1 Mpc -1 ) in all models. We used the luminosity 
function of Marzke et al. (1998) both considering a sin- 
gle Schechter function for all morphological types and 
three different Schechter fits for elliptical, spiral and ir- 
regular galaxies. Our counts are best fitted by considering 
three LFs, luminosity evolution and a moderate merging 
(Figure [l7| ) consistent with our interpretations. The merg- 
ing rate in Gardner (1998) follows Rocca-Volmerange & 
Guiderdoni (1990) with number evolution parameterized 
as (j>* oc (1 + z) v in the Schechter fit to the LF, and in 
order to conserve the luminosity density, L* cx (1 + z)~ n ; 
r\ is a free parameter. Good fits to our data are obtained 
with r\ < 0.5 in the case of three LFs for the cosmological 
parameters considered. 
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Fig. 17. Number counts models compared to various 
counts obtained in literature and in this work. All mod- 
els adopt Marzke et al. (1998) luminosity function, the 
dot-dashed line refers to three Schechter fits for elliptical, 
spiral and irregular galaxies with no evolution, the solid 
line uses the same LFs with luminosity evolution and a 
moderate merging, the dotted line is obtained by using 
a single Schechter fit without evolution, the short dashed 
considers luminosity evolution and the long dashed both 
luminosity evolution and strong merging for the same LF. 



8. Summary and Conclusions 

The HDF-S represents a unique opportunity for the study 
of faint galaxies up to now, both for its depth 
(F8UWabm™ ~ 29 for detection and F8WW A B,lim ~ 27 
for completeness) and spatial resolution (FWHMw0.2 arc- 
sec) . We presented here colors and number counts of HDF- 
S galaxies, along with number counts determined by split- 
ting the sample considering the morphology and the colors 
of the galaxies. We also analyzed the photometric data to 
constrain the redshift of HDF-S galaxies, and determine 
the contribution of different redshift populations to the 
counts. The main results are the following: 

— the number-counts relation has an increasing slope up 
to the limits of the survey in all four bands. The slope 
is steeper at shortest wavelengths; 

— the number counts model which best fits our data 
is obtained by considering three different Schechter 
functions for elliptical, spiral and irregular galaxies 
(Marzke et al. 1998), luminosity evolution and a mod- 
erate merging (77 < 0.5); 

— optical colors show that the sample contains a high 
fraction of galaxies bluer than local sources, for in- 
stance at F814Wab > 27 about 50% of sources have 
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(F450W— F606W)ab bluer than a typical local irreg- 
ular galaxy; 

— after an initial blueing trend the color-magnitude re- 
lation gets flat, suggesting a color-saturation due to a 
strong contribution of a flat spectrum population in 
the faintest bins; 

— morphological number counts ( F814WAB < 25) are 
dominated by late type galaxies = 0.43 ± 0.03, 
Jspir = 0.37 ± 0.03), while early type galaxies show a 
steep decrease in the faintest bin; 

— photometric redshifts of our sample galaxies show that 
the galaxies contributing with a steep slope to the 
number counts have z p hot ^ 1, suggesting a moder- 
ate merging; 

— morphological rcdshift distributions and number 
counts are biased against elliptical galaxies when using 
an apparent magnitude cut-off. 

Our number counts are in good agreement with previous 
results, except for the F300W-band number counts in the 
HDF-N. 

We do not see a flattening in F300W-band number 
counts, limiting our analysis at F300WytB <26.5. A steep 
F300W-slope may be explained both with a high fraction 
of low redshift galaxies or with merging or a mixture of 
them. In the former case, intrinsically faint, low redshift 
galaxies, imply a steep faint end slope of the LF in the 
F300W-band (a w —2) in the latter case we emphasize 
that if merging is the reason of the steep F300W-band 
counts, it would take place at z < 2 since the Lyman 
break prevents us from detecting galaxies in the F300W- 
band at higher redshift. A way of detecting merging and 
hence to check the validity of the hierarchical model, could 
be in principle the study of the evolution with look-back 
time of some structural parameters of galaxies, such as 
their intrinsic size. 

Morphological number counts are in agreement with 
previous results, as well as the morphological redshift dis- 
tributions. The morphological classification is, however, 
possible only for bright galaxies. The cut in apparent mag- 
nitude biases the sample against early-type galaxies, due 
to their large K — corrections. A correct test implying 
morphological classification should be done using K-band 
data (an extension of the test proposed by Kauffmann & 
Chariot, 1998) or a volume-limited sample, with morphol- 
ogy determined from spectra. Extracting a volume- limited 
sample from a magnitude limited survey does not fit the 
goal since at z < 1 almost 80% of the galaxies are too faint 
for a reliable morphological classification from imaging. 
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